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We propose a method for detecting the presence of synchronization of self- sustained oscillator 
by external driving with linearly varying frequency. The method is based on a continuous wavelet 
transform of the signals of self-sustained oscillator and external force and allows one to distinguish 
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of the signals. We apply the method to driven van der Pol oscillator and to experimental data of 
human heart rate variability and respiration. 
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I. INTRODUCTION 

It is well known that interaction between nonlinear 
oscillatory systems including the ones demonstrating 
chaotic behavior can result in their synchronization. Var- 
ious types of synchronization between oscillatory pro- 
cesses have been intensively studied in many physical, 
chemical, and biological systems pi 0, 0, 0, M, 0- Of 
particular interest in recent years is the investigation of 
synchronization in living organisms whose activity is de- 
termined by interaction of a great number of complex 
rhythmic processes. The underlying sources of these os- 
cillatory processes often cannot be measured separately, 
but only superpositions of their signals are accessible. 
For example, in electroencephalography recordings on 
the scalp the measured signals are the superpositions of 
the signals generated by various interacting sources. As a 
result, one can detect spurious synchronization between 
brain sources and come to wrong biological conclusions 
. Such situation is typical for many multichannel mea- 
suring devices. One faces the similar problem studying 
synchronization between the rhythms of the cardiovascu- 
lar system (CVS). The most significant oscillating pro- 
cesses governing the cardiovascular dynamics, namely, 
the main heart rhythm, respiration, and the process of 
blood pressure slow regulation with the fundamental fre- 
quency close to 0.1 Hz, appear in various signals: electro- 
cardiogram (ECG), blood pressure, blood flow, and heart 
rate variability (HRV) |8|]. This fact impedes studying 
their synchronization. 

Synchronization between the main rhythmic processes 
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in the human CVS has been reported in Refs. 0, 0, El 
[Til Il2| . It has been found that the systems generating 
the main heart rhythm and the rhythm of slow regulation 
of blood pressure can be treated as self-sustained oscilla- 
tors and the respiration can be regarded as an external 
forcing of these systems 0, However, at respira- 

tion frequencies close to 0.1 Hz it becomes difficult to 
distinguish the case of true synchronization between the 
respiration and the process of blood pressure regulation 
from the case of spurious synchronization caused by the 
presence of the respiratory component in the HRV and 
blood pressure signals used for the analysis of the rhythm 
with the basic frequency of about 0.1 Hz. Actually, the 
presence of external forcing can result in linear mixing 
of the driving signal and the signal of the self-sustained 
oscillator without any synchronization. Another possible 
case is the simultaneous presence of mixing of the signals 
and their synchronization. 

In this paper we propose a method for detecting syn- 
chronization of self- sustained oscillator by external driv- 
ing with linearly varying frequency. The method is based 
on a wavelet transform of both the external signal and 
the self-sustained oscillator signal and allows one to dis- 
tinguish the true synchronization from the spurious one 
caused by linear mixing of the signals. The principal 
moment of our approach is that we vary the frequency 
of external signal. We verify our method applying it to 
driven van der Pol oscillator and to experimental data 
of human HRV and respiration. It should be noted that 
synchronization between the different rhythmic processes 
can be detected also from the analysis of univariate data 
[ToL PHI fpH ITEl liH Il7j | . However, in the present paper we 
use only bivariate signals for detecting synchronization. 

The paper is organized as follows. In Sec. II the models 
for studying the effects of synchronization and mixing 
are considered. In Sec. Ill we describe the method for 
detecting synchronization using the continuous wavelet 
transform. Section IV presents the results of method 
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application to driven asymmetric van der Pol oscillator. 
In Sec. V the synchronization between the respiration 
and the process of slow regulation of blood pressure is 
studied. In Sec. VI we summarize our results. 



II. MODELS 

The universality of the phenomena observed in peri- 
odically driven self-sustained oscillators of physical and 
physiological nature has been discussed in Refs. [M Ho||. 
It has been shown there that qualitatively the same fea- 
tures of synchronization are observed in the case of peri- 
odic driving of van der Pol oscillator and in the case of 
respiratory forcing of the heartbeat and the process with 
the basic frequency of about 0.1 Hz. Let us consider the 
asymmetric van der Pol oscillator under external forcing 
with linearly increasing frequency as a model for study- 
ing interaction between the respiration and the process 
of blood pressure slow regulation. 

The periodically driven asymmetric van der Pol oscil- 
lator is described by the following equation: 

x — (j, (l — ax — x 2 ) x + fl 2 x — K sin <p(t), (1) 

where fi = 1 is the parameter of nonlinearity, £1 = 0.24-7T 
is the natural frequency, and K and tp(t) are, respectively, 
the amplitude and phase of the external force. The phase 

tp(t) = 2tt [(a + bt/T)] t (2) 

defines the linear dependence of the driving frequency 
u>d(t) on time: 

w«i(t) = ^ = 27r(o + 26t/T), (3) 
at 

where a = 0.03, b = 0.17, and T — 1800 is the maximal 
time of computation. We choose these parameter values 
to compare the results of simulation with those obtained 
at investigation of experimental signals of respiration and 
HRV (see Sec.©. 

The case a = in Eq. (JTJ corresponds to the clas- 
sical van der Pol oscillator with symmetric limit cycle. 
As a result of phase portrait symmetry, the power spec- 
trum of oscillations has only odd harmonics (2n + l)/o, 
n = 1,2,..., of the basic frequency fa. Since the sec- 
ond harmonic 2/o of the process with the basic frequency 
close to 0.1 Hz is well pronounced in the power spectrum, 
we consider the modified van der Pol oscillator with 
a = 1. 

The difference between the frequency luq of self- 
sustained oscillations and the natural frequency f2 is due 
to the effect of nonlinearity (/o = ujq/2h = 0.106 in 
Fig. (IJj). In the case of asymmetric van der Pol oscil- 
lator this difference is greater (/ = 0.098 in Fig.©). 

To compare the case of synchronization of oscillations 
by external driving with the case of mixing of the signals 
we consider the superposition of the signals 
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FIG. 1: (Color online) Power spectra and phase portraits for 
van der Pol oscillator Q at K — 0, f2 = 0.247T, and /i = 1 (a) 
a = 0, (b) a = 1 

where x(t) is the signal of the autonomous asymmetric 
van der Pol oscillator and R sin <p(t) is the additive sig- 
nal with the amplitude R, the phase ip(t), and varying 
frequency (J3J). 

III. METHOD OF DETECTING 
SYNCHRONIZATION USING CONTINUOUS 
WAVELET TRANSFORM. MEASURE OF 
SYNCHRONIZATION OF OSCILLATIONS 

Studying synchronization of chaotic oscillators vari- 
ous definitions of synchronization have been introduced, 
namely, the complete synchronization, generalized syn- 
chronization, lag synchronization, and phase synchro- 
nization Q. To investigate phase synchronization one 
has to choose the way of phase definition for the chaotic 
signals 0, III M M 

We will use the recently proposed approach |22ll23ll24l 
I25I ] to the analysis of synchronization based on examina- 
tion of a continuous set of phases defined with the help 
of the continuous wavelet transform |l8t |26J 

/+00 
x(t)r s , ta (t)dt (5) 
-00 



ais(t) = x(t) + i? sin 



(4) 



of the signal x(t), where ip s j (t) is the wavelet function 
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related to the mother- wavelet i(>o(t) as 



1 , (t-ta 
s V s 



(6) 



The time scale s corresponds to the width of the wavelet 
function ip Sj t (t) , to is the shift of the wavelet along the 
time axis, and symbol "*" denotes complex conjugation. 
We use the Morlet- wave let 271 



Mv) = (1/^tt) exp(jcrr?) exp (-?7 2 / 2 ) 



(7) 



as the mother- wavelet function. The choice of the wavelet 
parameter a — 2tt provides the relation s = 1// between 
the time scale s of the wavelet transform and the fre- 
quency / of the Fourier transform. 
The wavelet spectrum 



W(s,to) = \W(8,to)\exp\j<l> a (to)] 



(8) 



describes the system dynamics for every time scale s 
at any time moment to- The value of |W(s,to)| de- 
termines the presence and intensity of time scale s at 
the moment of time to- At the same time, the phase 
4> s (t) — arg W(s,t) is naturally defined for every time 
scale s. In other words, it is possible to describe the 
behavior of every time scale s using its phase 4> s (t) . 

If in the signals xi^(t) there is a range of time scales 
S\ < s < S2 for which the phase locking condition 



\A(f> s {t)\ = \<f> sl (t) - (f> s2 (t)\ < const 



(9) 



is satisfied and the part of the wavelet spectrum energy 
in this range does not vanish 
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(E(s)) ds > 0, 



(10) 



where (E(s)} is the distribution of integral energy by time 
scales defined as (E(s)) = (1/T) J t t+T \W(s, t )\ 2 dt , 
then the time scales s£ [s±; S2] are synchronized and the 
oscillators arc in the regime of time scale synchroniza- 
tion [23. In Eq. 10 (f> s i,2{t) are the continuous phases of 
the first and the second oscillator corresponding to the 
synchronized time scales s € [si;s2]- 

Using continuous set of time scales s and the phases 
associated with these scales we introduce the quantitative 
measure of chaotic synchronization p2l |25| : 



(E(s)) ds / / (E(s)) ds. 



(11) 



This measure defines the part of the wavelet spectrum 
energy fallen into the synchronized time scales. Increase 
of 7 from to 1 points to the increase of the part of the 
wavelet spectrum energy fallen into the synchronous time 
scales s. 
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FIG. 2: (Color online) Shaded plot of the wavelet power spec- 
tra |W(s,to)| for the external signal with the frequency vary- 
ing in accordance with Eq. J3J ( a )> the signal generated by 
oscillator Q (6), and the superimposed signal I^J (c). Time 
is shown on the abscissa and time scale is shown on the or- 
dinate. The color intensity is proportional to the absolute 
value of the wavelet transform coefficients. The scaies from 
the right side of the figure show the values of the coefficients. 
The dashed lines indicate the time scale Sd(t) corresponding 
to the linearly increasing frequency u>d(t) ©• 



IV. INVESTIGATION OF DRIVEN 
ASYMMETRIC VAN DER POL OSCILLATOR 

A. Amplitude dynamics of wavelet spectra of 
driven oscillator and superimposed signal 

Let us consider the wavelet power spectra |W(a,t)| of 
the external signal K simp(t) with linearly varying fre- 
quency, the signal x(t) generated by Eq. QJ, and the 
superimposed signal x^(t) Q. Typical wavelet power 
spectra of these signals are presented in Fig. [21 

The analysis of the wavelet power spectrum of the sig- 
nal x(t) of driven van der Pol oscillator (Fig. EJd) reveals 
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the classical picture of oscillator frequency locking by the 
external driving. As a result of this locking, the breaks 
appear close to the time moments t s and t2 S denoted 
by arrows, when the driving frequency is close to the 
frequency of autonomous oscillator or to its second har- 
monic. These breaks represent the entrainmcnt of oscil- 
lator frequency and its harmonic by external driving. In 
the region of the frequency entrainment the amplitude 
of the respective coefficients of the wavelet spectrum in- 
creases. This fact agrees well with the known effect of 
the oscillation amplitude increase in the synchronization 
(Arnold) tongue. If the detuning 8 — (toa — 27r/ ) is 
great enough, the frequency of oscillations returns to the 
natural frequency of the autonomous oscillator. 

It should be noted that besides the break at the main 
time scale so of the spectrum, the break at the scale sq/2 
corresponding to the second harmonic 2/o takes place. 
In this region the wavelet surface has no maxima asso- 
ciated with the driving signal since the intensity of the 
corresponding spectral component is low. 

The wavelet power spectrum of the superimposed sig- 
nal Xs(t) is shown in Fig. 0:. In contrast to the previ- 
ous case, the dynamics of both the van der Pol oscillator 
and the external force with varying frequency is well pro- 
nounced in the spectrum. The absence of break in the 
neighborhood of the time moment t s indicates the ab- 
sence of entrainment of the oscillator frequency. One can 
see only slight distortion of the wavelet surface caused by 
the increase of the wavelet coefficients amplitude \W\ at 
the time moments close to t s . This effect is the result of 
the addition of two signals with comparable amplitudes 
and close frequencies. Furthermore, the surface distor- 
tion in the region of the main scale sq is not followed 
by any change in the wavelet spectrum at the scale sq/2 
corresponding to the second harmonic of the signal x(t). 
Similarly, the coincidence of the frequencies aid and 47r/o 
at the time moment t2 S is not followed by any changes of 
the dynamics at the main scale so- 

Thus, it is possible to distinguish the cases of synchro- 
nization and mixing of the signals analyzing the dynamics 
of time scales corresponding to the basic frequency and 
its harmonics in the wavelet power spectrum. In the case 
of mixing the changes in the dynamics of the scale hav- 
ing the frequency which is close to the driving frequency, 
do not lead to any change in the dynamics of other time 
scales. In the case of synchronization the typical break 
in the wavelet power spectrum is observed at all charac- 
teristic time scales. 



B. Phase dynamics of driven oscillator and 
superimposed signal 

Let us consider the phase difference between the x(t) 
signal produced by Eq. and the external signal with 
linearly increasing frequency and the phase difference be- 
tween the Xs(t) signal produced by Eq. (0} and the ex- 
ternal signal. The phase differences A<p s (t) are calcu- 
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FIG. 3: (Color online) Phase differences between the super- 
imposed signal £Q an d the external signal with linearly in- 
creasing frequency (curves f and 2) and between the signal of 
driven van der Pol oscillator Q and the driving signal (curve 
S). R = 0.2 for the curve 1, R = 1.0 for the curve 2, and 
K — 0.2 for the curve S. The phase differences are calcu- 
lated at the scale Sd = 2n /uid- The fragment marked in (a) is 
enlarged in (b) 



lated along the scale Sd(t), corresponding to the linearly 
increasing frequency J3J), i.e., along the dashed lines in 
Fig. |21 Typical dependencies A(j) s (t) are presented in 
Fig. 

At first, we consider the case of linear mixing of the 
signals and analyze the phase dynamics in the vicinity 
of t = t s , where the driving frequency is close to the 
basic frequency of van der Pol oscillator: u>d(t s ) ~ 27r/o. 
The phase dynamics at the main scale sq = I /fa of van 
der Pol oscillator can be written as <j> s o(t) = "2irfot + 
(fo, where ipo is the initial phase. The phase <j) s {t) of 
the driving signal is determined by Eq. (J2J. Under the 
assumption that the amplitude of the external signal is 
significantly smaller than the amplitude of oscillations of 
van der Pol oscillator, the temporal behavior of phase 
difference between the components of the superimposed 
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signal xs(t) takes the form 

A</> a (t) = <pso(t) - Mt) = 2tt[(/o - a)* - (b/T)t 2 )} + (f . 

(12) 

It follows from Eq. I|12f) that the phase difference varies 
under the parabolic law and the parabola has extremum 
at t = t 8 . In the neighborhood of this moment of time the 
dependence A<fi s (t) is symmetric with respect to t = t s . 
Similar situation takes place at the scale sq/2 — l/(2/o) 
corresponding to the second harmonic. 

Such behavior of Acf> s (t) is illustrated by curves 1 and 
2 in Fig. [5] The curve 1 is plotted for small amplitude of 
the external force and has a parabolic shape in the neigh- 
borhood of t s and t2s ■ The curve 2 plotted for a greater 
amplitude R also has an extremum at t = t s . However, 
this extremum is less pronounced and has a shape differ- 
ent from quadratic one. It is explained by a great influ- 
ence of the external signal having the amplitude compa- 
rable with the amplitude of the signal generated by the 
van der Pol oscillator. This case is characterized by al- 
most constant phase difference which is close to zero at 
t far from t s . For large R values the part of the exter- 
nal signal in the mixture Xs(t) is sufficiently large. As 
a result, the phase difference is calculated between the 
phases of practically the same signals. Therefore, if the 
part of the driving signal is large enough in the super- 
imposed signal, one can detect spurious synchronization 
using conventional methods of synchronization investiga- 
tion. 

Let us consider now the case of synchronization of 
asymmetric van der Pol oscillator by external driving 
with linearly increasing frequency. This case is illustrated 
by the curve S in Fig. 

We analyze the behavior of the phase difference A(f> s (t) 
using the method of slow varying amplitudes under the 
assumption that the external driving does not change sig- 
nificantly the amplitude of nonautonomous oscillations 
but mainly has the influence on the phase relationships. 
It can be shown that A(f> s (t) is governed by Adler equa- 
tion m 



d(AMt)) 
dt 



n sin A(f> s (t) 



(13) 



-(27r/o -««,(*)) =0, 



where k is the coefficient defined by the oscillator pa- 
rameters. It follows from the Adler equation that in 
the region of synchronization defined by the condition 
(27r/o — u>d{t)) < K the phase difference A(j) s (t) monoton- 
ically decreases by ir under the driving frequency varia- 
tion. By this is meant that in the case of synchronization 
of oscillator by external driving with linearly varying fre- 
quency we will observe in the vicinity of time moments t s 
and t2 S the regions where the phase difference varies over 
7T. Such situation is illustrated by the curve S in Fig. 
A special feature of our approach for detecting syn- 
chronization is the application of the external driving 
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FIG. 4: (Color online) (a) Phase differences A<f> a (t) between 
the driven asymmetric van der Pol oscillator and the driving 
signal at scales Sd(t) + As. (b) Measure of synchronization 7 



with varying frequency. Another important feature of 
our method is the analysis of phase difference between 
the signals at different time scales s [2^, |2^| . We con- 
sider the behavior of phase difference A<f> s (t) at the scales 
s(t) = Sd(t) + As, where As is the detuning of the scale 
with respect to the basic scale Sd(t). 

In Fig. Uk, the phase differences at various scales are 
presented for the case when loq is close to UI4. It can 
be seen from the figure that for As £ (—1,2) the phase 
dynamics is qualitatively similar to the case of accurate 
adjustment to the basic scale Sd(i). At greater As val- 
ues the duration of epochs of synchronous behavior be- 
comes shorter and beginning from some values of detun- 
ing there is no synchronization at all. The presence of 
a range of scales As within which the synchronous dy- 
namics is observed allows one to investigate the system 
behavior without accurate adjustment of the scale to the 
basic scale Sd{t) varying in time. It can be useful for the 
analysis of experimental data. 

Figure^J) illustrates the dependence of measure of syn- 
chronization 7 defined by Eq. (|llfl on time. At t > 300 
the system quickly comes to the regime of time scales syn- 
chronization. The part of the wavelet spectrum energy 
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FIG. 5: (Color online) (a) Wavelet power spectrum to) | 

for the superposition of the signal of driven oscillator Q with 
K — 0.125 and the external forcing with R = 1. (b) Phase 
difference at the scale Sd(t) between the external forcing with 
R — 1 and the superposition of the external forcing and the 
driven oscillator Q signal with different K values 



fallen into the synchronized time scales is more than 80%. 
With further increasing of time the value of 7 decreases 
and at t > 540 it becomes equal to zero. The latter case 
corresponds to the absence of synchronization. 



C. Dynamics of nonautonomous system at 
simultaneous presence of mixing and synchronization 

Let us consider now the situation when the effects 
of mixing and synchronization of the signals of self- 
sustained oscillator and external driving are simultane- 
ously present. Figure shows the amplitude spec- 
trum |W(s,to)| of the wavelet transform for the signal 
being the superposition of the driven oscillator sig- 
nal synchronized by external driving with the amplitude 
K = 0.125 and the driving signal itself with R = 1. This 
spectrum has the features inherent to both the cases of 
synchronization and mixing and gives no definite answer 
to the question about the nature of the underlying dy- 
namics. The analysis of phase dynamics appears to be 
more informative. 

Figure \^jp shows the dynamics of phase differences 



Acj) s (t) for R — 1 and various amplitudes K. At small K 
values the situation is qualitatively similar to the consid- 
ered above case of mixing without synchronization. The 
regime of synchronization appears in a certain range of 
frequencies as the amplitude K increases. The variation 
of phase difference in the synchronization tongue takes 
the form of an inclined line. The boundaries of this line 
correspond to the loss of synchronization. Note that vari- 
ation of phase difference in the synchronization region in 
Fig. [SJs is smaller than in the case of R = 0. 



V. INVESTIGATION OF SYNCHRONIZATION 
BETWEEN THE RESPIRATION AND THE 
PROCESS OF BLOOD PRESSURE SLOW 
REGULATION 

This section contains the results of physiological data 
analysis. We studied eight healthy young male subjects 
having average levels of physical activity. The signals of 
ECG and respiration were simultaneously recorded in the 
sitting position with the sampling frequency 250 Hz and 
16-bit resolution. The experiments were carried out un- 
der paced respiration with the breathing frequency lin- 
early increasing from 0.05 Hz to 0.3 Hz. The rate of 
breathing was set by sound pulses. The duration of ex- 
periments was 30 min. 

Figure ®l shows a typical respiratory signal with lin- 
early increasing frequency and its wavelet power spec- 
trum (Fig. [U>). Extracting from the ECG signals the 
sequence of R-R intervals, i.e., the series of the time in- 
tervals between the two successive R-peaks, we obtain 
the information about the heart rate variability. Typ- 
ical sequence of R-R intervals for breathing at linearly 
increasing frequency is shown in Fig. [SJ:. Since the se- 
quence of R-R intervals is not equidistant, we developed 
a technique for applying continuous wavelet transform to 
nonequidistant data. Figure [SJi shows the wavelet spec- 
trum |iy(s,to)| for the sequence of R-R intervals pre- 
sented in Fig. |SJ:. This wavelet spectrum demonstrates 
the high-amplitude component denoted as 1 correspond- 
ing to the varying respiratory frequency manifesting itself 
in the HRV data. The second harmonic of the respira- 
tion is observed at a twice-higher frequency. The power 
of the rhythm with the basic frequency of about 0.1 Hz 
is rather small and this frequency is not pronounced in 
the spectrum. The amplitude of the respiratory rhythm 
in R-R intervals is much higher than the amplitude of 
the rhythm with the frequency 0.1 Hz. Comparing the 
wavelet spectrum of R-R intervals under linearly increas- 
ing frequency of respiration with the wavelet spectra pre- 
sented in Sees. Ill and IV one can come to a conclusion 
that a significant mixing of the considered physiological 
signals takes place without synchronization. However, 
the investigations performed in Ref. [Tl| with the same 
experimental data have clearly shown the presence of 1:1 
synchronization between the process of blood pressure 
slow regulation and respiration at breathing frequencies 
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close to 0.1 Hz. This synchronization has been observed 
within the time interval 200-600 s for each of the subject 
studied. The wavelet spectrum in Fig. EJl does not allow 
to detect the presence of synchronization. The increase 
of HRV amplitude at frequencies close to 0.1 Hz can be 
regarded only as an indirect indication of the presence of 
synchronization. 

In Fig. |7Jl the phase difference between the R-R inter- 
vals and respiration is presented. Between 200 and 600 s 
the phase varies in average almost linearly at the scale 
s = Sb(t) = l//b(t), where fb is the linearly increasing 
frequency of breathing. In this time interval the phase 
variation is close to tt indicating the presence of synchro- 
nization. Outside of the synchronization region the phase 
difference fluctuates around a constant value. From the 
result obtained it may be concluded that the respiratory 
dynamics manifested in R-R intervals affects the rhythm 
of blood pressure regulation within the time interval 200- 
600 s. Outside of this interval the respiratory component 
is observed in R-R intervals as a result of mixing of the 
signals without their interaction. 

Similar behavior of phase difference is observed at 
neighbor time scales. Figure [7jp shows phase differ- 
ences at the scales s = Sb(t) + As. For small detun- 
ing (As = ±0.75) the phase dynamics is qualitatively 
similar to the dynamics at the scale s — Sb- Hence, a 
small error in determining the required scale will not re- 
sult in qualitative changes in the plot. If the mismatch 
of the scale of observation and the scale of respiration 
Sb is large, the amplitude of the corresponding rhythm 
decreases and its phase cannot be defined confidently. It 
is the case As = —2 in Fig. 03. As a result, the phase 
difference is no more constant. 



VI. CONCLUSION 

We have shown that applying the external driving with 
varying frequency to the self-sustained oscillator and us- 
ing the method based on the wavelet transform one can 
distinguish the case of synchronization of oscillator by 
this driving signal from the case of absence of synchro- 
nization. The method allows one to detect the presence of 
synchronization between the signals with close frequen- 
cies even in the case when the effect of mixing of the 
signals is present. The use of the driving signal with vary- 
ing frequency is crucial for the proposed method, since 
calculation of the phase difference between the superim- 
posed signals at a fixed time scale gives no answer to 
the question, whether the system response is the result 
of active interaction of the oscillating processes or the 
result of their mixing without changes of the underlying 
dynamics. The proposed method does not require the ac- 
curate adjustment of the scale of observation to the time 
scale associated with the varying frequency of the ex- 
ternal driving. The efficiency of the continuous wavelet 
transform for estimation of the signal phase is demon- 
strated with experimental physiological data character- 
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FIG. 6: (Color online) Typical time series of breathing with 
linearly increasing frequency (a) and its wavelet power spec- 
trum (b). The respiratory signal is in arbitrary units. Se- 
quence of R-R intervals for the case of respiration with lin- 
early increasing frequency (c) and its wavelet power spectrum 
(d). The dashed lines indicate the time scale corresponding 
to linearly increasing frequency of respiration 
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FIG. 7: (Color online) Phase differences between R-R inter- 
vals and respiratory signal with linearly increasing frequency. 
The phase differences are calculated at the linearly varying 
scale Sb(t) = !//(,(£) (a) and scales si,(t) + As (b) 



ized by high level of noise. 

The proposed technique is applicable to the situation 
with superimposed signals being at the disposal, and it 
is necessary to define whether system 1 with adjustable 
frequency affects system 2 and can synchronize it. Such 
information can be useful for example for the diagnos- 
tics of the cardiovascular system state. Synchronization 
between the rhythmic processes in the human cardiovas- 
cular system under paced respiration is less effective in 
patients than in healthy subjects, and this effect corre- 
lates with the seriousness of the heart failure (2jJ. Using 
a small variation of the relative phase or high values of 
phase synchronization index as criteria of synchroniza- 
tion one can come to wrong biological conclusions due to 
the mixing of analyzed signals. 

The proposed method can be used for detecting syn- 
chronization in different systems if it is possible to vary 
the frequency of the external driving in the experiment. 
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